STATIONARITY OF PURE DELAY SYSTEMS AND QUEUES 
WITH IMPATIENT CUSTOMERS VIA STOCHASTIC 
RECURSIONS 



P. MOYAL 

Abstract. In this paper we solve a particular stochastic recursion in the sta- 
tionary ergodic framework, and propose some applications of this result to the 
study of regenerativity (that is, finiteness of busy cycles) and stationarity of 
some queueing systems: pure delay systems, in which all customers are imme- 
diately served, and queues with impatient customers. In this latter case under 
the FIFO discipline, we prove as well the existence of a stationary workload 
on an enlarged probability space. 



1. Introduction 

Real-time applications play a key role in the development of recent computer 
and telecommunications networks. In such a case the system must be able to 
meet severe delay requirements in the data transmission. The multimedia, time 
sensitive and on-line audio/video data traffic flows are subject to such temporal 
constraints, as well as call centers in which customers hang up if not attended after 
a short waiting time. Then it may be better to discard the requests which would 
eventually arrive too late, in order to favor some other packets which still can meet 
their delay requirements: they can therefore be thought as "impatient" since they 
may leave the system forever if not treated fast enough. 

These systems are commonly referred to as real-time queueing systems, and rep- 
resented in a random environment by queues with impatient customers: this is the 
classical stochastic model of a queueing system fed by a stationary process of ar- 
rivals, marked by the service times of the customers, enriched by a new (possibly 
random) parameter, the initial patience or time credit of the customers. Therefore 
the customers are lost forever when their patience expires before they could reach 
the service booth (in this case the patience is until the beginning of service), or 
even complete their service (it is then said to be until the end of service). The time 
credit of the customer decreases at unit rate as time goes on, and provided that the 
customer is in the system at t, his residual time credit at t is the remaining time 
before his deadline, i.e. the time already spent in the system substracted to the 
initial time credit. 

A first question arising in the study of random dynamical systems, and in particu- 
lar of queueing systems, is the one of stochastic stability. In the general stationary 
ergodic framework (i.e. for G/G queues), in most cases this problem is handled 
by studying at arrival or departure times some stochastically recursive sequence 
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{F n } neN , that is, such that for some initial conditions, 

Y n +1 = 4>n(Yn), 

where <j> n is a stationary ergodic sequence of random maps. In "good cases" (es- 
sentially, when the almost sure continuity and monotonicity of <j> n are assumed), 
backwards recurrence schemes as in |Loy62| exhibit a stationary version of {Y n } nefsS , 
which can be finite and unique depending on some condition on the parameters. 
In the latter case, provided that the stationary solution is possibly nonpositive and 
that the system is empty on {Y n — 0}, the system is regenerative, i.e. it empties 
almost surely an infinity of times. 

It is a straightforward fact, that under general G/G/ assumptions the queue 
with impatient customers is almost surely finite in infinite time, since its conges- 
tion process (i.e. the process counting the number of customers in the system) is 
dominated in a suitable sense by that of a well chosen pure delay system: there are 
an infinity of servers, so that every customer is immediately served upon arrival. 
The size of this latter system is known to be a.s. finite in the long run (see (BB02J, 
Property 2.8.1.), which implies in turns the same property for the queue with im- 
patient customers. But regenerativity or existence/uniqueness of some stationary 
state are all but obvious properties even for trivially finite systems (see in [FH83 
counter-examples for the workload sequence of the G/G/ 1/1 queue, which may 
admit either none or more than one stationary version), and at the time, no such 
strong results as have been proven for both queues with impatient customers and 
pure delay systems. 

We believe that this lack of information is due to the strong instability of the 
system's dynamics. Consider for instance a queue with impatient customers treated 
under some service discipline depending on the time credits of the customers, say 
Earliest Deadline First (EDF for short), in which case the servers chose when 
available the customer with the smallest residual time credit. This discipline has 
been shown to be optimal for feasability (see |Der74| l. instataneous loss for strong 
stochastic orderings ( |PT88] ) and loss probability in the Poissonian case ( [t>M03 
Moy05))). In such a case, the situation (waiting time, vocation to be lost or served) 
of an arriving customer not only depend on the past of the system, but as well 
on the subsequent arrivals of prioritary customers. Therefore, so far no exhaustive 
description by a simple stochastically recursive sequence has been presented in the 
litterature. Remark, that representations by counting measure-valued processes 
have been proposed to study the mean flow of this queue (see |DM05| . DM06 
or |DLSf)l| l. where Dirac measures are affected to each residual time credit at 
current time. 

Hereafter we present original descriptions of pure delay systems and queues with 
impatient customers at arrival times by remarkable real- valued sequences, all of 
them satisfying a particular recursive equation (equation (E) below). The study 
of the stationary regimes of these sequences then provide sufficient information to 
yield conditions of regenerativity and convergence to a stationary state for these 
systems. 

To go one step further, in the particular case of queues with impatiences treated 
under the FIFO (first in, first out) discipline, we state that the workload sequence 
(waiting time proposed to the arriving customers) is stochastically recursive (see 
equation ltl9|l below), but in this case the recursive functional (j> n is not monotonic. 
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We circumvent this difficulty using the recent techniques of enlargement of proba- 
bility space (see |AK97j ). which lead to the existence of a stationary version of the 
sequence in a weaker sense. 

This paper is organized as follows. In section [21 we investigate in the stationary 
ergodic framework the existence and uniqueness of a solution for the stochastic 
recursion (E). We present applications of this result in section to the study of 
stationarity of pure delay systems, in section^lto queues with impatient customers 
until the beginning of service, and in section|^lto queues with impatient customers 
until the end of service. In subsection 14.21 we show the existence of a stationary 
workload in the case of a single-server proceeding under the FIFO discipline. 



2. Stationary solution of a stochastic recursion 

Consider a probability space (il, T, P°) , embedded with the measurable bijective 
flow 9 (denote 6> -1 , its measurable inverse). Suppose that P° is stationary and 
ergodic under 9, i.e. for all 21 £ T, P° [tf- 1 ^] = P° [21] and (921 = 21 implies 
P° [21] — or 1. Note that according to these axioms, all ^-contracting event 
(such that P° [2l c n 9~ X A\ = 0) is 0-invariant. We denote for all n £ N, 6 n = 
9 o 9 o ... o 9 and 9~ n = 9~ 1 o 9~ l o ... o 9~ x , and say that two sequences of r.v. 
couple when there exists a P°-a.s. finite rank N such that they coincide for all 
n > N. 

Let a > 0, f3 > be two nonnegative integrable random variables (r.v. for 
short), and the sequences of r.v. [a„, /3n}„ e z = {ao6 n ,/3 o 9 n } neZ . Let Z be an 
a.s. finite r.v.. Consider the following stochastically recursive sequence : 

= Z, 



(E) 



yZ 

yZ 
I n+1 



= [max {Y£ , «„} - fin] + for all neN. 



In this section we are interested in the existence and uniqueness of a stationary 
solution {l^} ngN to the recursion (E), in that we are looking for an a.s. finite r.v. 
Y such that for all n € N : 

Y 7 Y = Yo9 n . 
In particular, Y would be a solution to the equation 
(1) Yo9= [max{r,a}-/3] + . 

We have the following: 

Theorem 1. There exists an only P°-a.s. finite solution Y of {^P - given by : 



(2) 



Y := 



sup 



j 

E 



Moreover, for all P°-a.s. finite and nonnegative r.v. Z, the sequence {^^} neN 
couples with {Y o 9 n } n£N . 

Proof. Existence We exhibit a solution Y of with a backward recurrence scheme, 
following the ideas of the construction in |Loy62j. Define the following sequence: 

Mo = 0, 



max a_ 



j 

E 

i=l 



for all neN. 
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Remark, that {M„} n£N is non decreasing, since for all n S N : 



Mt+i = max < M„, a_„+i - ^P-i? 



Thus, the P°-almost sure limit of {M„} ngN is given by : 



On another hand, for all n € 



sup a - 3 



i=l 



(3) M n+1 o . 





f » / 


max < 


max 







, a_i - /3_i 



i=0 



iax a_j - /3_ t 



max 

3 



o9-\ a_i l-/3_i 



= [max{M ra , a} — (3] + . 
The function [max{., a} — (3] + o _1 being P°-a.s. continuous, we thus have : 

Moo o0= [max {Afoo, «} - /3] + , 

hence Moo is a solution of (QJ. Now, Birkhoff ergodic Theorem ( [RobOO], Theorem 
10.9) entails that: 



3 



-E° [/?] , P°-a.s. 



Thus for all e > 0, there exists N e 6 N such that : 



P u [6J := P 



> 1 -e. 



i=l 

For e < min {E° [/?] , l}, on £ e , for all j > AT e , 

a_,-^/?_ 4 <j(-E° \J3]+e) < 0, 

i=l 

hence Afoo < oo. Therefore : 

P° [Moo < oo] > 1 - e > 0. 
But as easily seen, {Moo = oo} is 0— invariant, hence : 

P° [Moo < oo] = 1. 

Therefore, Moo provides an a.s. finite solution of 

Uniqueness Let Z be an P°-a.s. finite solution of JQ. On the one hand, Z > = Mo 
P°-a.s. and Z > M n P°-a.s. implies: 

Z o = [max {Z, a}~ I3] + > [max {M n , a} - /3] + = M n+1 o , P° - a.s., 
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and thus Z > M n+1 P°-a.s. .Therefore, for all n £ N, P° [Z > M n ] = 1 and to the 
limit : 

(4) P [Z>M oo ] = l. 
On another hand, we have : 

(5) P° [Z < a] > 0. 

Indeed, if P° [Z = 0] > 0, there is nothing to prove. Now if P° [Z = 0] = 0, 
P° [Z < a] = entails that on a P°-almost sure event, 

Zo9>0o [max {Z,a} ~ (3} + > 4^ Z o9 = max {Z, a} - /3 = Z - j3. 

We thus have : 

E° [Z o 9 - Z] = -E° [f3] < 0, 

which is absurde in view of the ergodic lemma ( |BB02j . Lemma 2.2.1.): J^J is 
verified. Since 

{Z < a} = \Z o 9 = [a - /?]+} C {Z o 9 < o 9} , 
© implies that 

(6) P° [Z < Moo] = P° [Z o 9 < Moo ° 6] > 0. 
But on {Z < Moo}, 

Zo9= [max {Z, a} - j3} + < [max {Moo, a} - (3} + = Moo ° 0, 
which means that the event {Z < Moo} is ^-contracting. Thus, with Jfj} : 

P° \Z < Moo] = 1, 

and with <@J : 

P° [Z = M oa ] = 1. 
Let us give the developped form of Y® for all n € N : 



(7) 



Y° = 



max | a,- 

j'=o 



n-l 



with the conventions rnax^Q . = and Ylj=a ■ = ®- Indeed, F ° = and assuming 
Q for some neN leads to: 



x n+l — 



max ■ 



n-l 

max | a,- 

3=0 



Ea 



max < max I a 7 - — 3i \ ,a n — B n 

. 1 J ' \ I 

Let < Z < oo, P°-a.s. .Remark, that for all n e N, 



max | a 

3=0 



*=3 



(8) 



i=0 
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as easily seen by induction: 



Y Z = Z = max{Z,0}, 



and if (jHJ holds for some n £ N, 



i=0 



n-1 
, max I a 




= max <^ - ]T 

Thus, 

{y„ z ^ y„° for all n € N} = {Y z > F n ° for all n £ N} 

= <^ r r f = Z - Pi > for all n e N 
But the latter event is of probability zero since it contradicts 

-E° W] < 1 , 



( lim -fz-Vft 

I \ i=0 



which is P°-almost sure from Birkhoff 's theorem. Hence, since two sequences driven 
by the same recursive equation couple from the instant when they coincide for the 
first time, there exists a random indice 

t z <° = inf {k e N, = Y° for all n > k} < oo, P°-a.s. 

This is true for any finite and nonnegative r.v. Z, and in particular for Z = 
Y. Thus, for all rank larger than t z > y := max {t z >° , t yo } < oo, {Y z } and 
{Y Y — Yo 6 n } neN coincide: these two sequences couple. □ 

Corollary 1. For all finite and nonnegative r.v. Z, there exists P°-a.s. an infinity 
of indices such that Y„ = if and only if : 



(9) 



sup a-j - /3-i < 



i=i 



>0, 



Proof, in view of the previous coupling property, 

n r t 

i—l i—1 i—1 

= P° [Y = 0] , P°-a.s.. 



Hence 



E 1 {y-=o} = °°. P °- a - s - P ° i Y = 0] > 0. 



which is condition Q. Now if © does not hold, on a P°-almost sure event, 
the sequence j^rf } ngN coincides after a certain finite rank with the R+*-valued 
sequence {Y o #™} ngN , and hence never reaches after this instant. This event is 
9ti -invariant, and thus P-almost sure. □ 
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3. STATIONARITY of the G/G/oo system 

We hereafter consider a pure delay G/G/oo system : under general stationary 
ergodic conditions, customers enter in a system where there is an infinite number 
of servers. They are consequently immediately served upon arrival, and spend in 
the system a time equal to the duration of their service. It has been known for 
a while (see for instance BB02J, p. 133), that this system is stable, in the sense 
that its congestion process (number of customers in the system) is a.s. finite after 
an infinite time. In this section, we address the question of stability in terms of 
existence and uniqueness of a stationary version for some process describing the 
queueing system (here the LRST process, see below). This study will as well lead 
to a condition for the system to be regenerative, i.e., to empty a.s. an infinity of 
times. 

On the probability space (f2, !F, P), we therefore define the stationary ergodic 
point process N, whose points {T n } ngZ represent the arrivals of the customers 
{Cn}„ e z> with the convention that To is the last arrival before time t = 0. We 
denote (#t) igR , the associated measurable bijective flow, such that for all t € M and 
all 03 e B(R): 

N(0 t u)(*8) = + t), P-a.s., 

where 23 + t := {x + t,x e 23} . The interarrivals, denoted for all n £ Z by £„ = 
T„ + i — T n , hence form a stationary sequence of integrable r.v., and we denote 
£ := Co = Ti — Tq. The process N is marked by the stationary sequence {cr n } n£Z 
denoting the service durations requested by the customers, and the integrable r.v. 
cr := (To is the service duration of customer Co. 

Let us denote (O, T ", P°, 0), the Palm probability space of N(a), where 9 := #Ti is 
the associated bijective flow (see for example |BB£)2|). In particular, P° is stationary 
and ergodic under 8, T = P°-a.s. and for all n e Z, 

£,n=£o9 n and a n = aod n . 

Since there is an infinite reservoir of server, the customers are immediately at- 
tended, and for instance C n spends in the system the time a n necessary for his 
service. Provided that T n < t < T n + a n (i.e. C n in the system at time t), the 
remaining processing time (that is, time before the end of service) of C„ at this 
instant is given by T n + <j n — t. Let us denote for all t € R, 

%t — ^ 1 {T„<t<T„+cr„}, 

the number of customers in the system at t. Provided that %t > 0, we denote 

m 1 (t)<m 2 (t) < ...<m Xt (t), 

the remaining processing times of the customers in the system at t, ranked in 
increasing order, and 

2) t :=m Xt (t)Mx t >o} 

the largest remaining processing time (LRPT for short) of a customer in the system 
at t (zero if the system is empty). 

Let us depict the evolution of the LRPT process. If the system is empty just 
before the arrival time T n (2)t„- = 0), then 2)r„ = c„, i.e. the amount of service 
customer C n requires. If not (2)t„- ^ 0), then 2)t„ is the maximum between 
2)t„- (remaining service time of the customer who will stay in the system the 
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latest among those present at T n —), and the remaining service time of C„ at T n , 
i.e. cr„, Inbetween instant of arrivals, (2)t) t>0 decreases at unit rate: at a departure 
time (end of a service), either there are still customers in the system, in which case 
the value of the LRPT is striclty positive, and doesn't jump at this time, or the 
customer leaving the system was alone in it, i.e. the LRPT equals zero at this 
instant. This dynamics is thus described by the following equation: 

(10) 2) t = [max{2) T „_, a n } - (t - T n )} + , t G [T n ,T n+1 ] . 

Note, that (2)t)teK has right continuous with left-hand limits (rcll for short) paths, 
and that 11 Oil implies for all n G Z: 

2)t„ + i- = [max{2) Tn _,a„} -£„)]+, 

which is a stochastic recursion of the type of (E) on (ft,P°). We aim to find a 
(possibly unique) r.v. 2)o such that for all t£t: 

(11) 2) f ( w )=2)o(M,P -a-s.. 

The process (2)t)t>o wom d in particular be a stationary process, being compatible 
with (0t) t>0 - We also address the question of the regenerativity of the queue, i.e., the 
almost sure infinity of instants in which the queue empties. We have the following: 

Theorem 2. There exists an only P°-a.s. finite solution 2)o of given by : 



(12) 2)o := 

Moreover, provided that 



sup (cr„ + T n ) 

n<0 



(13) 



sup (a n +T n ) < 

n<0 



>o, 




the system is regenerative. 

Proof. From Theorem there exists an unique finite r.v. 

Y := 
such that: 

(14) Y o0 = [max{y,cr} -^] + , P°-a.s.. 

We thus can construct the stationary LRPT sequence {^T} neZ an d the corre- 
sponding stationary LRPT process (2)t) teK , defining for all neZ: 

=Yoe Tn , 

*Q t = [max {Y^,a n } - (t - T n )] + , t e [T n , T n+1 ) . 
From CH, the event 



lim 2) t for all n 6 Z 



jF^j = [max {F/, cr„} - $n] + for all n G z} 
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is P°-almost sure and ^-invariant, hence P°-almost sure. Thus, P°-a.s. (2)t) tgK is 
rcll, and satisfies iflOj) . 
Now, since 



2} = [max{Fo^ To ,(7} + To]" 
for all n G Z, t S [T n , T n+ i), 



i + 



sup (a„ + T n ) 

n<0 



P°-a.s.., 



2)o ° Ot = [max {F o To o t , a o 6 t } + T Q o t ] + 

= [max{2) T „_, <r„} - (t - = 2) t , 

the process (2)t) tgR hence satisfies (11 II . On the other hand, denote for all non- 
negative and finite r.v. Z, {Y^ } neN '■= {?)f„-}„ e N ! LRPT sequence initially 
equal to Z. Provided that ltl3|l holds, we can apply Corollary^ 

P° [there are an infinity of indices n such that Y% = 0] =1, 
which implies, since the latter event is 9t ± -invariant, that 

P [there are an infinity of indices n such that 2)f„- = 0] =1. 
The system is thus regenerative since for all t, 

{X t - 0} = {2} t - 0} . 

□ 



4. Queues with impatient customers until the beginning of service 

Following Barrer's notation ( |Bar57j ). we now consider a queue with impatient 
customers until the beginning of service G/G/s/s+G(b). Customers enter in a 
system with s non-idling servers and of infinite capacity, according to a stationary 
ergodic point process (N t ) t>0 of points {T„} ngN , and as above we denote {£n}„ e N 
the stationary sequence of interarrivals. The customers request service of durations 
{c n } n6N , and are furthermore impatient until the beginning of their service, in that 
they leave the system if they do not reach the service booth before a given deadline : 
customer C n is discarded from the system if the time he must wait in the buffer is 
larger than D n (his initial patience, or time credit). He is thus eliminated forever, 
at time T n + D n . The sequences {<7n} neN and {D n } n£N are sequences of integrable 
and non-negative marks of {N t ) t>0 . We denote again for all t > 0, Xt the number of 
customers in the system at t. The servers follow a non-preemptive service discipline: 
we therefore consider that a customer can not be eliminated anymore as soon as he 
enters the service booth, even though he meets his deadline during his service. 

We define the Palm probability space (fi, T , P°, 9) of the arrival process N(a, D), 
using the same notations as above. 



4.1. Regenerativity in the general case. 
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. According to the asumptions made above, the total sojourn time of customer 
C n does not exceed D n + a n , i.e. the sum of his total patience and the time 
necessary for his service. On the other hand, it is at least equal to min {<j n , D n }, 
i.e. the time needed for him to be lost, or immediately served. Hence, provided 
that C n entered the system before t (T n < t) and even though he already left 
the system before t, his remaining maximal sojourn time at t (i.e. the remaining 
time before his latest possible departure time, if not already reached) is given by 
[er„ + D n — (t — T n )] + , whereas his remaining minimal sojourn time at t (i.e. the 
remaining time before his earliest possible departure time, if not already reached) is 
given by [min {a n , D n } — (t — T n )] + . Hence the largest remaining maximal sojourn 
time (LRMST for short) at t among all the customers entered before t is given by 

C t := max [er„ + D n - (t - T n )} + 

n—l 

and the largest remaining minimal sojourn time (LRmST for short) at t, by 
l t := max [min {er„, D n } - (t - T n )\ . 

n—l 

These two processes evolve according to the following dynamics: they decrease at 
unit rate inbetween arrival times (there is no possible change among customers 
already entered), and equal the initial maximal (resp. minimal) sojourn time of C n 
at his arrival time T n , provided that it is larger than the value of the process just 
before T n . In other words, for all n S Z and all t 6 [T„, T n +i): 

(15) C t = [max{£ T „-,(T„ + D n } - (t - T n )] + , 



(16) t t = [max{[ Tn _,min{er„,L>„}} - (t - T n 

Theorem 3. The G/G/s/s+G(b) queue is regenerative if: 



i + 



(17) 

and only if: 
(18) 



sup (a n +D n +T n ) < 

n<0 



>o, 



sup (min {a n , D n } + T n ) < 

n<0 



> 0. 



Proof. Define for all finite nonnegative r.v. Z and all n 6 N: 

L z n :=C Tn - 



l z 



the LRMST (resp. LRmST) just before the arrival of customer C„, provided that 
£t - = Z (resp. It - = Z). The processes (£t) tgK and (k) teR have rcll paths, 
hence in view of equations l(T5|l and ifTfijt. we have the following recursive equations: 

L n+1 = [ max { L fn °n + D n } - £„] + , 
l n+l = [ max { l ni min Wn,D n }} - £„] + . 

On (il, P°) the two latter equations are of type (E), hence similarly to the argument 
of the previous proof, corollary implies that there are P-a.s. an infinity of indices 
such that if = if and only if itTTjl holds, and an infinity of indices such that i„ = 
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if and only if ltl8|l holds. We conclude remarking that for all initial conditions and 
all t > 0: 

{A = 0}C{£ t = 0}C{[ 4 =0}. 

□ 

4.2. The G/G/l/l+G(b)-FIFO case. 

. Let us now consider the special case, where there is one server proceeding ac- 
cording to the FIFO (First in, first out) order. Denote for all t 6 1, 2Ut the 
workload submitted to the server at time t, i.e. the quantity of work he still has 
to achieve at this time, in time unit. The process (2U t ) tgR (already investigated 
in JBH81 and [BBH84 for the renewal case) has rcll paths, and we define for all 
n, W n = 2U T -. Its value at t equals the work brought by the customers arrived 
up to t, and who will eventually be served, since the other ones won't ever reach 
the server. Under the FIFO discipline, the served customers are those who find a 
workload less than their patience upon arrival. Inbetween arrival times, (2U t ) teR 
decreases at unit rate. Hence, for all n G Z and t G [T„, T„+i): 

W t = [2U T „_ + cr n l {mSTn _< Dn} - (t - T„)] + , 
whereas the workload sequence obeys the following recursive equation 

(19) W n +i = [W n + a n l {Wn < Dn} - £„] + . 

For all n G N, let W® be the workload seen by C„ upon arrival, provided that 
C enters an empty system (that is, W® = 0). We have the following: 

Lemma 1. For all n G N, 

W°<L° n , P°-a.s.. 

Proof. For all n G N, on the event {W° < 

• If D n < W n and L° < D n + cr„, then 

L n+1 = + D n - U + > [W° - £»] + = W° +1 . 

• IfD n > W° and < D n + a n , we have 

L n+1 = [a n + D n - £„] + > [W° n + a n - £„] + - W n ° +1 . 

• If D n < Wn and L° > D n + a n , then 

L n+1 =[L° n -^] + >[W n -^] + = W n+v 

• Finally, D n > W° and L n > D n + a n yields 

L°n+1 = K ~ Cn] + > [Dn + (T„ - £„] + > + (7„ - £„] + - TU° +1 . 

Therefore, in any case L n+1 > since Lg = Wq P°-a.s., by induction L„ > 

W n for all neN, P°-a.s.. □ 

On (f!, P J , let </){■} be the measurable random map from (R, Q5(M)) into itself 
(we denote 4> G 9Jt(R)) defined by 

{y} = [y + a(Lu)i {y < D{uj)} - £(^)] + . 

Equation iflflll can thus be rewritten as : 

(20) Wn+i =<f>(6 n .){W n }, P°-a.s.. 
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In addition to the previous result of regenerativity (Theorem |2J| , we investigate in 
this particular case the existence of a stationary solution for the stochastic recursion 
(1201 . i.e. a finite r.v. W such that W o 9 = <fi{W}. As opposed to equation 
(QJ for instance, the construction of such a solution is complexified by the non- 
monotonicity of the function <j> in the state variable, in which case a construction 
of Loynes's type becomes fruitless. Nevertheless, it turns out to be possible to find 
a "weaker" stationary solution, that is, defined on an enlarged probability space. 
To this end, we follow the steps of the construction presented in |A K97: : on the 
product probability space O x R, we define the shift 9 by 

9{u,x) = {9Lj l( j){uj){x}) . 

We then have: 

Theorem 4. The stochastic recursion admits a weak solution, that is, a 9- 
invariant probability P° on fl x R whose fl-marginal is P . Therefore, on (fi x R) 
there exists a R x 9Jt(R) -valued r.v. ( W, <j>) satisfying: 



Wo9 = 



In particular, is stationary under P , and has the 

same distribution as {</> o #™} ngN . 

Proof. From Theorem^! {L° } N couples with {L o 0"}„ eN , where 



L := 



SUp U-j 

jew V 



D. 



Therefore this sequence tends in distribution to Y: it is in particular tight, which 
implies with Lemma that {W r °} neN is tight, since for all e > 0, there exists M e 
such that for all neN, 

P° [W° < M e ) > P° [L° n < M e ] > 1 - e. 

Define now on ft x R the random variables : 



W(u, x) := x, <f> (uj, x) := <j)(w), 



and for all n G 



W n (uj,x) :=w(e n (w,x) 
Remark that for all n G N, 21 G T and 58 £ S(R), 

= P° ® 5 [21 x 



and 
P°®<5 



P°®c5o 

(n x 05) 



<T" (21 x IE 
= P° ® <J 



P°[2l] 



0-" (w^C®)) 



= P° <g> (So W - -^®) = P° [W° G 58] 



Hence, the probability distributions <^ P° <g) S a o 0- n \ on Q x R have f2-marginal 

I J neN 

P° and R-marginals the distributions of {Wn} ne n) which form a tight sequence : 
the sequence \ P° ® So ° 0~ n > is tight. We can therefore apply Theorem 1 

I J neN 
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of |AK97| : there exists a ^-invariant probability P° on x R whose Sl-marginal is 
P°. It is now straightforward that : 

W n {(J, x)=(j) (e^uj) o <f> (9 n - 2 uj) o ... o <j>{oj){x}, 

4>o8 n (u>,x) =(f>oO n {w), 
hence the sequence \ W n > satisfies on f2 x E the stochastic recursion: 

W n+1 = $o§ n {w n }, 

where \ W n , 4>oQ n \ is stationary under P°. □ 

I J n&l 

5. Queues with impatient customers until the end of service 

. Let us now consider a G/G/s/s+G(e) queue : the model is that of the previous 
section, except that the customers are now assumed to remain impatient until the 
end of their service. Indeed, they leave the system, and are eliminated forever 
provided that their service is not completed before their deadline. Keeping the no- 
tations and other assumptions of the previous section, customer C n is thus dicarded 
when the total time he has to wait in the buffer and spend in the service booth is 
larger than his initial time credit D n . To avoid the trivial case, where a customer 
can not even wait for the end of his own service (in which case he is automaticaly 
rejected), we furthermore assume that D n > a n , for all n G Z. 

In this easily seen the LRMST and LRmST processes are respectively 

given for all t by: 

C t := max^ [£>„ - (t - T„)] + 
[ t : = max^ r i 1 [a n - (t - T n )} + , 

and obey the following dynamics: for all n £ Z and all t € [T n , T n+ i) , 

C t = [max {C Tn _,D n }-(t-T n )} + 
l t = [max{[ Tii _,cr„} - (t - T n )] + . 

Hence, similar arguments as those of section ITTl lead to the following result: 
Theorem 5. The G/G/s/s+G(b) queue is regenerative if: 



sup (D n +T n )<0 

n<0 



>o, 



and only if: 

P° sup(<7 n +T n ) <0 >0. 

n<0 

On the other hand, in the special case G/G/1/1+G(e) under the FIFO discipline, 
a customer is served provided that the workload he finds upon arrival doesn't exceed 
the time he agrees to wait in the buffer, i.e. D n — a n . Hence, the process (2Ut) teR 
satisfies for all n G N and t G [T„, T n +i): 

W t = [W Tn - + "Vi 1{2IJt„ - <Z) n — CTn } - (t - T n )] + , 

and the workload sequence satisfies the recursive equation: 
(21) W n+ i =iP(6 n .){W n }, P°-a.s., 
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where 

V>(w) {y} =[y + ^)^{y<D(^)-a(^)} - . 

We define the following shift on SI x i: 

8(w,x) = {du,ip{u)){x}). 

As above, denote for all n G N, W^ 1 (resp. L°) the workload (resp. LRMST) 
found by customer C n provided that the system was empty at Tq — . Then, for all 
n W„ < lP n , as easily checked taking D n := D n — a n in Lemma Hence Theorem 
1 of [AK97| can still be applied and leads to: 

Theorem 6. The stochastic recursion 121]) admits a weak solution P° on fl x M. 
In particular, there exists a R x 9Jl(W) -valued r.v. (W,ip) satisfying: 

W o9 = ${W} , 

{(W, o 9 n } nGN is stationary under P°, and {tjj o #"} neN ftas ifte same distribu- 
tion as {V> O ^"Ingpj. 
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